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Abstract A theory is derived relating the configuration of the cores of active 
O i galaxies, specifically candidates for presumed super-massive black hole binaries 

<sJ '< (SMBHBs), to time-dependent changes in images of those galaxies. Three de- 

^ ' flection quantities, resulting from the monopole term, mass quadrupole term, and 

spin dipole term of the core, are examined. The resulting observational technique 
is applied to the galaxy 3C66B. This technique is found to under idealized circum- 
stances surpass the technique proposed by Jenet et al. in accuracy for constrain- 
^ ■ ing the mass of SMBHB candidates, but is exceeded in accuracy and precision 

^SJ I by Jenet's technique under currently-understood likely conditions. The technique 

OO . can also under favorable circumstances produce results measurable by currently- 

available astronomical interferometry such as very-long baseline-interferometry 
(VLBI). 
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1 Introduction 



The deflection of light by gravity is the oldest experimentally-verified test of the 
' theory of general relativity [T] . With the continued improvement in observational 

resolution in astronomy, particularly through very-long-baseline interferometry 
(VLBI), the detection of more subtle effects of this light deflection becomes prac- 
tical. Consequently, light deflection can be used to measure the properties of dis- 
tant systems. This paper supplies a theory for using time- variable light deflection 
to measure or constrain the parameters of binary systems. Specifically, the de- 
flection angle of a light ray from a distant source is related to the configuration 
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and motion of a binary system located in a distant galaxy somewhere between the 
point of emission of the light ray and its observation. 

Super-massive black hole binaries (SMBHBs) are thought to form the cores 
and primary energy sources of the broad class of galaxies termed "active galax- 
ies", "blazars", or "quasars". However, a combination of distance, radio noise, 
and optical thickness makes direct observation of presumed SMBHBs impracti- 
cal. Observing a time-dependent motion in the image of the galaxy can provide 
information on the mass and orbital parameters of an SMBHB candidate. 

Work by Damour and Esposito-Farese |6| and by Kopeikin et al. [5| estab- 
lishes a theory of time-dependent light deflection by describing the time- 
dependent part of the deflection through the quadrupole term, which is the lowest- 
order term resulting from the mass distrubtion whose effects are practical to eval- 
uate using current astronomical observational techniques. The work of Mashhoon 
and Kopeikin |25 | in examining gravitomagnetic effects furthermore provides a 
theory for evaluating the contribution of the spin dipole of such systems and com- 
plements the work of Einstein 1 27 1 in providing a complete theory for stating the 
location of the deflected image, in the weak field limit. We generalize these theo- 
ries to a stronger-field regime and put constraints on the theory's applicability in 
this regime. 

As a case study of an active galaxy, the theory is applied to the galaxy 3C66B, 
a nearby active galaxy with a candidate SMBHB core II13II . and theoretical con- 
straints on 3C66B's parameters from a light deflection experiment are compared 
to the constraints claimed by Jenet et al. |2j. 

2 Theory 

2. 1 Notations, definitions & assumptions 

We assume that Einstein's theory of general relativity is true to the limits of our 
ability to observe and applicable to the systems under examination. We do not 
address MOND or other post-Einsteinian models. 

Throughout this paper, "emitter" refers to the source of light rays being ob- 
served; "deflector" refers to the mass distribution causing a change in the metric 
of spacetime from flatness; and "observer" refers to the point where the light rays 
produced by the emitter are observed. 

We also make use of a coordinate system derived from the Cartesian system, 
defined thus: in a space that is asymptotically Cartesian let a line be described bjQ 

x\t)^k'(i-t^)^x\,. (1) 

Let t*\>t the time associated with the line's closest approach to the origin of the 
Cartesian system. Let t — t — t* denote a new time coordinate (that is, at r = 
the line reaches the closest point to the origin of both the Cartesian and projected 
systems). Space coordinates are projected onto a plane passing through the origin 
of the coordinate system and perpendicular to a line from the observer to the origin 
of the coordinate system; these new space coordinates are denoted ^* = P'^^x^ {t*) 

' Throughout this paper, Greek indices indicate (0,1,2,3) and Latin indices indicate 
(1, 2, 3). A contravariant 3-vector is denoted either in boldface or with a raised Latin index. 
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where the projection operator is defined P'^ = (5'^ — fc*fc^ . In the projected coor- 
dinate system, the index refers to r and the indices i denote coordinates 

For a trajectory described by ([T]| let = P*^ (t) be the "vector impact 
parameter" of the trajectory and let d = |^*| be the "scalar impact parameter" 
of the trajectory. Since the space is asymptotically flat, d is also the ratio of the 
magnitudes of the angular and linear momenta o f the light ray. Note then that 
for the trajectory described by (r), r (r) — \J (P ^ Let the unit vector 
= C7d- 

We assume that the wavelengths of all light rays observed are much shorter 
than the longest wavelength of gravitational radiation emitted by the deflecting 
system. 



2.2 General theory 
2.2.1 Background 

Consider a photon emitted at some distant point at some time in the distant 
past to- This beam of light in asymptotically flat space follows a path such 
that the coordinate a;* of the photon is given by the relation ([T]) ; therefore, /c* = 
{dx^ ld£)\^__^. Let fc* be normalized such that k^ki — 1; then the vector fc" = 

(l, fc') is parallel to the four- momentum of the photon in flat space. 

Let an asymptotically-flat metric be a function of some affine parameter 
A. Let " = fc" + (A) + S"" (A) be the trajectory of a photon moving in 
this metric space, where describes the part of the trajectory arising from the 
spherically-symmetric non-flat part of the metric and S"" describes the trajectory 
arising from a perturbation to the metric. Then, we have the geodesic equation 13] 
equation 87.3] 

''^^"^^^"^ +r^",i^^g^^o. (2) 

The quantity d (k" + E°-) /d\ corresponds to the change in momentum of the 
light ray in space, which when projected onto a plane of observation corresponds 
to the angular deflection of the light ray. We define this deflection vector by |j5l|6l 

"^(*'0^^'k■ + -^]observer + ^«^ (3) 

where the term Z\a' corresponds to corrections arising from any contribution to 
deflection other than our deflector. 

In the case of Eddington's experiment O on solar deflection, the "true" po- 
sition of the emitter - that is, the position of the emitter observed in the limit of 
intervening deflection going to zero - was known. In the case of deflectors with 
small proper motion, in this case extragalactic or otherwise distant objects, where 
the emitter would be seen without the presence of the intervening deflector may 

^ Throughout this paper we use the convention G = c = 1 to simplify our equations. 

' All metrics in this paper are stated using the harmonic gauge condition, that is, 
g^" r^v ~ 0. The Minkowski metric in Cartesian coordinates is chosen with signature 
(— , +, +, +) and is denoted r]^,j, and we make use of the Einstein summation convention. 
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not be known; therefore, the periastron of the light ray must be determined by 
other means. Let P be the periastron of the light ray's trajectory about the deflec- 
tor; in such cases, the time delay between the deflection and the motion of the 
deflector is related to the periastron by 

^ ~ ^peak deflection ^ '^alignment' 

where tpg^k deflection tims when the image of the source is observed to 

be deflected most from the position of the deflector and ^alignment time 
when the projected images of the components of the system and the image fall 
into a line, assuming that P < p/2 and that the change in the gravitational field 
propagates at the speed of light. 



2.2.2 Description of deflector 



Our deflector of interest is as follows: two objects are denoted with the indices 
1 and 2. The mass of object 1 mi > 7712. The objects have positions x\ (t) and 



(t) and velocities vl (t) and vl (t) [FIGURE [H. Then our source has density 



x: 



distribution 



p (t^ x^) = mi5 (x' - x\ (t)) + 17128 {x^ - x\ (t)) 
and velocity distribution 

(t, x^) ^ v\8 [x^ - x^^ + v\8 [x^ - x^^ 



(5) 



(6) 



where 5 (x*) is the three-dimensional Dirac delta distribution. 

Our metric has the form = 77^;^ + s^j, + h^^y where s^J_y is the non- 
Minkowski part of the Schwarzschild metric and h^^y is a small perturbation. Let 
hf^jj be the perturbation resulting from the quadrupole moment of the mass distri- 
bution and let h^^^ be the perturbation resulting from the spin dipole of the mass 
distrubtion. Let the variable s — t — r. Then explicitly, the metric is given by 16] 
|29l|25]|26lp 181]: 



(7) 





2m 
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Sij 






y,2 ' 



and 



'00 



dx^dx^ 



°'' "dx'dt 



(8) 



dx^dx^ r 



dt"^ r ' 
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and 



"■00 ^ "•ij ^ "i 



(9) 



where the vector S"' = ( J^^, J^i, Ji^) and J*^' = / (xT^" - x^T^o) 
chapter 2.9]. 

Let the objects orbit one another with a known period p. Let our coordinate 
system origin be located at the center of mass of the binary and let m = toi + m2. 
Let = x\ — x\ ^ vector denoting the spatial separation of the two masses and 
1= I a' I . Let the mass ratio q = ni2 / nii < 1. 

By our choice of coordinates, the dipole term of the deflector's mass distribu- 
tion is zero. 

Using the Landau-Lifschitz definition of the transverse traceless 
quadrupole |[3] eq 41.3], the quadrupole moment of the deflector is given by: 



Qij (0 = I P (x, t) [ZxiXj - r^5ij\ dV = — ™^ 2 [^'^'^ 
J (1 + 9 



I Si. 



(10) 



In the case that the masses travel in almost circular orbits about their common 
center of mass, then in a primed coordinate system related to our chosen system 
only by unitary rotations. 



a" {t) = I 



/sin(2^ + 0')\ 




5a" [t) 



(11) 



where 0' represents a constant phase term, and where Ja'* is small. Rotating from 
the primed system first about the y-axis, then the x-axis, then the z-axis, we have 



a' (t) = I 



( 



\ 



cos tf^ sin 



sin ^ sin 



sin If' sin 6? cos + (/)^ \ 
I- cos If' sin Q cos 



27rt 
V 



COS 6) COS ( ^ 

V 



5a' (t) 



J 



(12) 

where cf) subsumes rotation about the y-axis with (/>' and where and if are the 
angles of rotation of the plane of motion away from the jcz-plane about the x- and 
z-axes respectively. 



2.2.3 Solution to the geodesic equation 



The theory of the effects of small perturbations to the metric on light propagation 
in the weak-field limit is already developed |5 6|. However, since the effects of a 
quadrupolar perturbation fall off as d^, it is desirable to expand the theory to be 
applicable to regions of stronger fields. We note in particular that for a closely- 
orbiting compact binary system such as a SMBHB, then m and / will be of similar 



magnitude; therefore, we extend the first-order theory of light deflection to order 

Oim/df. 

First, note that all terms in (|7]l are O (m/r) or higher and that all terms in ([8| 
are of O (mP/r^) < O {rn? /r^). Let O [m? /r^) be small such that all higher 
orders are negligible. Then, suppressing negligible terms. 



hi- -2 



(13) 



Let the Christoffel symbol associated with the Schwarzschild metric 7^^^ ' = 
— (1/2) (77"'' + s"'') {sps,-^ + s^s.p — S/37,<5) and the remaining part resulting from 

the perturbation F^^^ = - (1/2) (??"'') {hf3s,^ + h-^s^p - hp^^s)- Then © be- 
comes 

r;f ^ + ^) {k^ + ^^ + s") {k-^ + + s^) = . (14) 

Since all F^}^ and all components of must be at least of O {m/r) or zero, 
( 114b expands, again suppressing negligible terms, to 

+ + r;f^ (k^ + {k^ + K?) + rf^h^k^ = . (15) 

Since 



r;^^^^(fc'5 + «'^)(A;^ + «;^) = 0, (16) 



we conclude 



_^^a^(h)^/3^7 ^ (17) 

which is exactly the result for the weak-field approximation |[5]|6l. 

Plugging (II6I 1 and (117b into ([31 and choosing r as our affine parameter, we can 
define the Schwarzschild and non-Schwarzschild(Q!jj) parts of the deflection 
angle[FlGURE|2l: 

a\, (D ^ P^H^ (18) 

and 

a\ {t, C) = P;S' = / ihps,j + h,sj3 - V,^) k^k-^dr. (19) 

The monopole term a\,i (^') of the deflection produced by the core is static and 
unique, regardless of changes of configuration within the core |i4jJJ. We can use 
the general, exact solution for provided by Darwin 11611 : 

For purposes of this derivation only, choose spherical coordinates. By the sym- 
metry of the monopole term, this part of the trajectory of the light ray must lie in a 
plane, so we can choose the coordinate 6 as an affine parameter and the coordinate 
<j) as constant. Then we obtain an equation of motion 
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Identifying the impact parameter with a conserved quantity in the system 
(r^/ (r — 2m)) [dt/d9) — d and substituting in, have three solutions; we discard 
the two where the light ray never reaches a distant observer and take the remaining 
one, 

1 Q-P + 2m Q -P + 6m 



r (0) 4mP AmP 



-sn^C {0) , (21) 



where the constant Q is defined by = {P — 2m) (P + 6m), the periastron and 
impact parameter are related by d^ = P^ / {P — 2m) and 

C {6) = y^iQ/P) (0 + 9o), and snC is the Jacobi elliptic sn function fBO] 16.1.5]. 
In the limit of P 3> m, inverting this relationship and taking its asymptotic limits 
at large r leads to the well-known relationship 

4:771 

"M,weakfield(eO = — (22) 
As P ^ 3m, however, the deflection becomes 



fj. (f ) = In ^ + In 648 (tVS - 12 



TTl 

In— + 1.248 (23) 

d 



where /i is the angle of deflection about the apse of the trajectory, rather than the 
deflection seen by a distant observer; the angles involved are no longer necessarily 
small so we cannot approximate aM = /i. In the case of an impact parameter 
comparable to 3^7, it is no longer useful to consider the monopolar displacement 
in and of itself as small differences in impact parameter cause great changes in 
deflection angle, and multiple images of a source may be detectable, some of 
which may result from geodesies which travel several times around the deflector. 
Our consideration therefore must focus not on the static deflection but on time- 
dependent deflections arising from higher multipole moments of the deflector. 

Kopeikin and Mashhoon |25 1 develop the effect of the rotation of a system on 
that system's deflection of light, in the weak field approximation. Investgation of 
this effect is useful for the system as described in that every practical case of an 
astronomical binary will display orbital motion. However, the theory developed 
by Kopeikin and Mashhoon is only sometimes compatible with the strong-field 
approximation presented herein. 

The integration of ^ is trivial. Let Let [C) be that part of determined 
by h^^ . when the deflector is stationary relative to the observer, the resulting con- 
tribution is given by 

«S (r) - ^ {2S^k''n'ejkin' + k'S'^e^.k) ■ (24) 

Calculating 5* with (II2I) for the case of a binary whose components are in 
almost-circular orbits, 

f, ,2 / 9 „\ /sin'f'cos©, , 

= ^m^— cos cose, . (25) 

1 + 9 P \ + {i+qf-{2.i/pf \ -sine 
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We must emphasize that (l24l) is compatible with the O (m'^/r^) generalization 
above only when O {mP /dPp) > O (m^/d'^); in particular, when Z « p the sys- 
tem's motion is no longer slow. We draw attention to this contribution to empha- 
size the difficulty in asssociating an image with a particular source and to under- 
score the utility of time-dependent deflection versus time-independent deflection 
in parameterizing a system. 

Let aQ^(t,i^*) be that part of a\ determined by h^jj. Uq is determined by 
plugging ® into ( 1191 ); while |6| uses the method of Fourier transforms, the form 
of ^ allows direct integration of a Fourier series decomposition as well; either 
way, the result is the following deflection |5 60: 

^0 = StT^ [("2 - (*)) - «i (•^) «2 (*) e^fefc^'^l (26) 
" (1 + 

for which we reiterate the following properties: firstly, in contrast to the monopolar 
case where aJw always points along <^', the quadrupolar deflection has a contribu- 
tion parallel to^% aQ\\, and also a contribution perpendicular to aq^, which 
vanishes only in the case that a component of the projected quadrupole moment 
vanishes, that is, only if the axis of rotation of the deflector is perpendicular to 
our line of sight; and secondly, the deflection depends only on the configuration of 
the deflector at the time of the light ray's closest approach to the center of mass, 
t = t*. 

In the case of almost-circular motion, inserting ( 1121) into ( 126b leads to 



i COS 2^ 



(l + sin^ e) cos(^^ + 2(t) 
+ sin^ 0~l 



- sin 2W sin O sin f ^ + 2^ 



and 



it,e 



mq 



isin2tf' 



<P (1 + qf 



(-i + sin^ e) 

+ (l + sin^ e) cos + 
cos 2^ sin 6* sin + 20^ 



2(j) 



> n" 



(27) 



e .,7c 



Pn't28) 



The relationships ( 127b and ( 1281) are original to this work and have not previ- 
ously appeared. From this relationship it is easy to see that the time-dependent 
deflection of the emitter's image is periodic, with a period half that of the orbit of 
the core's components. 

The greatest time-dependent deflection is observed when the emitter lies on 
the line of the semimajor axis of the apparent motion; when q — 1; and when the 



The symbol eijk represents the Levi-Civita permutation symbol defined such that ei23 = 1. 



9 



plane of the system lies perpendicular to the plane of observation. In this case, 
( |27] | reduces to 



3/2 

ag {t, d) < 



4:7rs 

cos ( h 20 ) - 1 



(29) 



so the total quadrupolar deflection seen over one half -period of the deflector's 
motion is 



3P 

Aaq (d) < --^m. 



(30) 



Compared to the monopole deflection in the case of a large impact parameter, 



OLM 



< 



3P_ 
Ad?' 



(31) 



2.2.4 Other contributions to the deflection angle 

If the path of the light ray after its closest approach to the deflector but far from 
the deflector is nearly occulted (for example the Sun or another star), then deflec- 
tion from this intermediate deflector, a? . ( f' • . |, must be taken into 

mtermediate v^njnty' 

account as well. Where^* refer to the vector impact parameters of light rela- 



tive to these intermediate deflectors, dn = C ■ ^ 
deflectors, and where m„/(i„ is small for all n. 



, and nin are the masses of these 



"intermediate v^nJ j j ■ y^'^' 

In linear approximation and in the harmonic gauge, the various deflections 
can be superposed linearly. The total deflection of the light ray from our source, 
therefore, is given by 

{t, e) = «b e) + "M ie) + a^s ie) + ^intermediate (C,int) • (^3) 



2.3 Application to 3C66B 

3C66B, also known as 0220h-43, is a radio galaxy l9\ with z = 0.0215 |10|, ap- 
proximately 91 Mpc distant from the Milky Wajjj. 3C66B exhibits jets emerging 
from its core, making it a good candidate for the location of a SMBHB 1 12 1. 



^ We use a value of 71 km/s/Mpc for the Hubble constant Ho for all distance calculations. 

im 
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2.3.1 Parameters of the system 

Sudou et al. |131 give estimates of 3C66B's parameters by direct radio observa- 
tion of its core, including an upper limit on m, a period, and an orientation of 
the core's motion. Sudou also reports a limit on the minimum impact parameter 
available for determining the parameters of the system using a first-order approxi- 
mation theory, corresponding to the limit of optical transparency at VLBI's higher 
operating frequency. The parameters Sudou gives are: 

m < 4.4 (1 + qf x 10^" solar mass = 6.5 (1 + qf x lO^^ cm, 
I < 5.1 (1 + g) X 10^^ cm, 
P > 23^ as = 3.1 X 10^^ cm, 
d > 3.7 X 10^^ cm, 
p= 1.05 ±0.03 years, 
O = 15° ±7°, 



(34) 



where P is constrained by the limit of the core's opacity in the radio spectrum and 
is derived from the apparent eccentricities of the elliptical boundaries of radio 
opacity. From I and P we can furthermore conclude that in the case of maximized 
l,q< 0.20. 

2.3.2 Estimates for distant emitters 

Although highly eccentric motion in 3C66B is not ruled out 11411 . the age of the 
presumed binary is great enough to have circularized the orbit through gravita- 
tional radiation under most conditions 1 15 1. We present the case of circular motion 
as an upper limit on the time-dependent deflection angle, noting that if all other 
parameters are constant then in the case of eccentric motion any time-dependent 
separation of the masses must have / as an upper bound. 

Using the maximal figure for mass and the minimal figure for impact parame- 
ter in ( 134b and applying ( 1231) , the ratio ni/l — 0.30, placing our proposed system 
in the regime of strong deflection. We find a monopolar deflection of: 

/ fi 5 X 1 0^^ cm \ 

u = In (1.44) + 1.248 = 0.13 radian = 7.2°. (35) 

' V3.7xl0i6cm^ 7 

The components of the system as proposed by Sudou have 27r//p < 0.39. 
Therefore it is not reasonable to apply ( 1241) to 3C66B in the regime where deflec- 
tion from the quadrupole moment will be detectable. 

Deflected images lying along the major axis of the core with the system as 
constrained in (134b will have time-dependent deflections in the following amounts: 

12(5.1x10^'' Cm)^(l. 2)^ (6.5x10^^^ Cm)(1.2)^(0.2) , , ^g) 

Z\q!qII (d) < (3.7x1016 cm)^ (Tip j . 

— 5.8 X 10~^ arcsecond 



AaQ^ (d)<i||!^(0.26) 

■ ^(1-2)^ (6.5xl0i^ cm 

1.4 X 10 ^ arcsecond 



12(5.1x1016 cm)^(i.2)^ (e.sxiQi^ cm)(i.2)^(o.2) (37) 
AaQ±_ (d) < (3.7x1016 cm)-' (Tip (0.2b) 
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with a period of p/2 — 0.53 ± 0.02 yearsfor each component of the deflection. 



3 Observational techniques 

3.1 Interferometry 

Electromagnetic interferometry provides the best currently-available techniques 
for high-resolution astronomy. The use of space-based interferometry and im- 
provements in equipment allowing for higher frequencies of observation continue 
to steadily improve resolution capabilities. The current most powerful technique 
available is VLBI, which Sudou et al. used to determine the motion in the core of 
3C66B |13|. 

VLBA, the Very Long Baseline Array, is an array of ten radio telescopes II19II 
operating in wavelengths as short as 3mm operating as a single large interferome- 
ter. The current best available resolution is 



1.7 X 10^^ arcsecond |f20|, making VLBA currently capable of constraining the 
parameters of 3C66B further through direct observation as well as the Jenet pul- 
sar timing experiment described below accomplishes indirectly. The launch of the 
space-based ASTRO-G satellite CIJ will extend the resolution capabilities further. 

The SIM PlanetQuest mission (formerly Space Interferometry Mission), cur- 
rently scheduled for launch in 2015 [22], is expected to have a resolution capabil- 
ity of 4 X 10^^ arcsecond [23 1. SIM will operate in the optical band and quasar 
observation is part of the planned mission. 

Farther into the future, the MAXIM (Micro-Arcsecond X-ray Interferometry 
Mission) satellite array currently in development 1 24 1 is expected to give resolu- 
tions on the order of 10~^ arcsecond in the x-ray band, and is explicitly designed 
with the observation of black holes in mind. 



3.2 Pulsar timing 

Jenet et al. 12 examined the period of the pulsar PSR B 1855-1-09 for changes in 
its period over several years, motivated by the idea that as gravitational waves 
generated by the core of 3C66B pass through the pulsar then the pulsar's signal 
should be modulated with a period related to the period of the proposed 3C66B 
SMBHB. The distance between the Solar System and the pulsar furthermore give 
the advantage that the signals observed modulating the pulsar are some 4000 years 
older than the motion observed in the 3C66B core, providing some information 
about the evolution of the system. However, Jenet' s experiment produced a null 
result. 

The experiment's analysis involved examining the frequency space of the pul- 
sar's signal for components in a range from (1/27.8) yr^^ to 19.5 yr^^, then sub- 
tracting out the one-year and six-month components resulting from geodetic ef- 
fects. The results are described as showing no signal distinguishable from noise 
other than the already-known main oscillation frequencies of the pulsar. Therefore 
the magnitude of gravitational waves generated by 3C66B, and consequently the 
parameters of its core, can be further constrained. 
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Jenet et al. assert that a system with 

m {q/ {I + q) \ > 0.7 X 10^° solar mass can be ruled out by the observed 

null resuh in the change in pulsar periods over seven years; this corresponds in 
the optimal case of g = 0.2 to a system with m = 2.3 x 10^° solar mass = 
3.4 X 10^^ cm. For a system under these new constraints, we estimate optimal 
peak deflections: 

12(5.1X10^° Cm)^1.2)'^ , 15 . (0.2) „„^ 

AuQw {d) < (3.1x101^ cmf 1^-4 X 10 cmj ^ (1.07j . 

= 2.1 X 10 ^ arcsecond 

12(5.1X10^^ Cm)^(1.2)^ . 15 , (0.2) . . 

Z\aQ_L (d) < (3.1x1016 cm)^ l"^-^ ^ (Tlf ^^'^^^ 

= 5.0 X 10~^ arcsecond 

which remains within the detection limit of VLBA as currently configured. 



4 Conclusions 

A theory of light deflection by time-dependent distributions of matter has been 
presented for metrics which are perturbations of the Schwarzschild metric, ac- 
counting for deflection resulting from time-independent and time-dependent terms 
in the metric. To order m^/r^, deflections originating from the quadrupole mo- 
ment of the mass distribution and, with some constraints, the dipole moment of 
the system's spin can be linearly superposed on the system as if in a weak-field 
approximation. The theory can be practically evaluated for and applied to a model 
of the core of an active galaxy, but the theory of light defiection from the spin of 
the deflector needs further development for applicability in the regime of strong 
deflection. 

The examination of time-dependent light deflection is a feasible technique for 
the evaluation of proposed SMBHB systems, imder idealized circumstances. In 
the event that a suitable emitter exists, examination of light deflection can be used 
to constrain the parameters of the proposed SMBHB in the core of 3C66B. We 
emphasize that while the existence of an identifiable suitable emitter in the case 
of 3C66B is unlikely, the theory can be applied equally well to any other SMBHB 
candidate, any of which may have a suitable source; in particular, active galax- 
ies displaying Einstein rings or other artifacts of strong gravitational lensing, es- 
pecially multiple images, should be examined. The theory can be equally well 
applied to intragalactic objects. 

The quadrupolar motion in the core of 3C66B can be examined and evaluated 
by the observation of deflected images in the region of the sky near the core of the 
galaxy, if found; the time-dependent part of the deflection has a magnitude of up 
to 58 microarcseconds parallel to the impact parameter of the emitter and up to 14 
microarcseconds perpendicular to the impact parameter. 

For the case of 3C66B, for most emitters pulsar timing can constrain the pa- 
rameters of the deflecting system better than time-dependent light deflection can. 
VLBA in its current configuration is capable of constraining the parameters of the 
core of 3C66B under ideal circumstances. Anticipated interferometers will have 
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resolutions up to two orders of magnitude greater and will be capable of evaluat- 
ing the parameters of the system more closely while examining it in a wide range 
of frequencies, and may make the observation of time-dependent light deflection 
resulting from motion in the core of 3C66B more practical. 
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S' 




Fig. 1 Object 1 has mass mi , velocity vl and is located as position x\ ; object 2 has correspond- 
ing m2, V2 and X2. x\ — X2 = a* and the spin vector 5* where S' = (j^^, J^"', J^^^ and 
J'^ = / [x^T^° - x^T^°) dV is perpendicular to a\ vl and v^. 




Fig. 2 A light ray produced by the emitter initially follows trajectory A:\ which has its closest 
approach to the origin of the coordinate system at In a pure Schwarzschild space, the light 
ray follows trajectory fc' + (A) and is detlected about the apse of its trajectory by angle //; in 
a perturbed Schwarzschild space, it follows trajectory fe' + (A) + (A) and the observer 
records an additional deflection of a\. 
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